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Abstract 



Let A = A* he a linear operator in a Hilbert space H. Assume that equation Au = 
f (1) is solvable, not necessarily uniquely, and y is its minimal-norm solution. Assume 
that problem (1) is ill-posed. Let fs, ||/ — < S, be noisy data, which are given, 
while / is not known. Variational regularization of problem (1) leads to an equation 

I A*Au + au = A*fs. Operation count for solving this equation is much higher, than for 
solving the equation {A + ia)u = fs (2). The first result is the theorem which says that 
if a = a{S), lim5_^o ^(5) = and lim^^o = 0, then the unique solution us to equation 
(2), with a = a{6), has the property lim^^o I |'W<5~1/| I = 0. The second result is an iterative 
method for stable calculation of the values of unbounded operator on elements given with 

'. an error. 

I 1. Introduction 

The results of this note are formulated as Theorems 1 and 2 and proved in Sections 
1 and 2 respectively. For the notions, related to ill-posed problems, one may consult PP 
and p] and the literature cited there. 

Let A = A* he a. linear operator in a Hilbert space H. Assume that equation 



is solvable, not necessarily uniquely, and y is its minimal-norm solution, y-LN := {u : 
Au = 0}. Assume that problem (1) is ill-posed. In this case small perturbations of / may 
cause large perturbations of the solution to (1) or may throw / out of the range of A. 
Let fs, 11/ — < S, he noisy data, which are given, while / is not known. Variational 
regularization of problem (1) leads to an equation A*Au + au = A*fs, where A* = A 
since we assume A to be selfadjoint. Operation count for solving this equation is much 
higher, than for solving the equation 



Au = f, 



(1) 



{A + ia)u = fs 



(2) 
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The result of this paper is the following theorem. 

Theorem 1. Let A = A* be a linear bounded, or densely defined, unbounded, self- 
adjoint operator in a Hilbert space. Assume that a = a{6) > 0, lim^^o ci,{^) = and 
lim5_^o = 0, then the unique solution us to equation (2) with a = a{S) has the prop- 
erty 

lim 1 — ?/| I = 0. (3) 



Why should one be interested in the above theorem? The answer is: because the 
solution to equation (2) requires less operations than the solution of the equation {A*A-\- 
al)u = A*fs basic for the variational regularization method for stable solution of equation 
(1). Here / is the identity operator. Also, a discretized version of (2) leads to matrices 
whose condition number is of the order of square root of the condition number of the 
matrix corresponding to the operator A* A + a{6)I, where / is the identity operator. 

Proof of Theorem 1. One has 

\\ua,s-y\\ < \\{A + ia)-\fs-f)\\ + \\{A + za)-'Ay-y\\ < ^ + a\\{A + ia)-'y\\. (4) 



a 



Moreover, 



\ima^\\{A + ta)-'Ay-y\\^ = ]ima'' T ^^^^ = 0, (5) 

a^O a^O J_^ + 



where we have used the spectral theorem, Es is the resolution of the identity correspond- 
ing to the selfadjoint operator A, and we have taken into account that 



lim a' / "r^"'"^ = 



° d{Esy,y) 



a^O 7_g + 

because y-LN. From formulas (4) and (5) one concludes that if a = a{S) > 0, lim^^o = 
and lim^^o = 0, then the unique solution us to equation (2) with a = a{S) satisfies 
equation (3). 

Theorem 1 is proved. □ 
2. Calculation of values of unbounded operators 

Assume that A is a densely defined closed linear operator in H. We do not assume in 
this Section that A is selfadjoint. If / G D{A), then we want to compute Af given noisy 
data fs, 11/5 — /II < S. Note that fs may not belong to D{A). The problem of stable 
calculation of Af given the data {fs,S, A} is ill-posed. It was studied in the literature 
(see, e.g., by a variational regularization method. Our aim is to reduce this problem 
to a standard equation with a selfadjoint bounded operator < B < I, and solve this 
equation stably by an iterative method. 

Let V = Af. This relation is equivalent to 

Bv = Ff, (6) 
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where B := {I + Q)~^, F := BA, Q = AA* is a densely defined, non-negative, selfadjoint 
operator, the range of / + Q is the whole space H, and S is a selfadjoint operator, 
< B < I, where the inequalities are understood in the sense of quadratic forms, e.g., 
B>0 means {Bg, g) > for all g e H, and F := (/ + Q)-^A. 

Lemma 1. (see 0) The operator {I + Q)~^A, originally defined on D{A), is closable. 
Its closure is a bounded, defined on all of H linear operator with the norm < ^. One has 
{I + Q)~^A = A{I + T)~^ , where T = A* A is a non-negative, densely defined selfadjoint 
operator, and ||A(/ + r)~^|| < i. 

If fs is given in place of /, then we stably solve equation (6) for v using the following 
iterative process: 

= (/ - B)v^ + Ffs, Vo±N\ (7) 

where A^* := {u : A*u = 0}. Let y be the unique minimal-norm solution to equation (6), 
By = Ff. Note that y = Hy + Ff, where H := I - B. 

Theorem 2. If n = n{6) is an integer, lim5_^o '^('^) = oo and lims-,o[Sn{5)] = 0, then 

lira - y\ \ = 0, (8) 

where vs := Vn{s), and f„. is defined in (7). 

Proof of Theorem 2. From (7) one gets Vn+i = J2^=o H^F fs + H'^'^^uq, where 
H := I — B . One has y = Hy + Ff. Let Wn '■= Vn — y- Then Wn = Yl]Zo H^Fgs + H^Wq, 
where g^ := fs — f, and Wq is an arbitrary element such that Wq-LN*. Since < H < 
-^^ll-^ll < |, and \ \gs\ \ < one gets 



Wr. 



<'^ + [j\l-sf-d{E,,wo.w,)Y'\ (9) 



where Eg is the resolution of the identity corresponding to the selfadjoint operator B. 
If WqLN\ then 



1 

2n 



lim / {l-sy''d{Es,Wo,w^) = \\Pwo\\\ (10) 



h 



where P is the orthoprojector onto the subspace {u : Bu = u} = {u : Qu = 0} = N* , 
and Pwq = because wo-i-N* by the assumption. The conclusion of Theorem 2 can 
now be derived. Given an arbitrary small e > 0, find h sufficiently close to 1 such that 
(1 — s)^^d{Es,wo,wo) < e. Fix this h and find n = n{6), sufficiently large, so that 
6n{5) < e and, at the same time, (1 — hY'^^^^ < e. This is possible if 6 is sufficiently 
small, because lim5_^o ^('^) = oo and lim5_^o[^'^('^)] = 0. Then (1 — sY^d{Es,wo,wo) < 
JI^{1 — sy^d{Es, wq, Wq) + J^{1 — sY''^d{Es, wo, Wo) < e, and inequality (9) shows that (8) 
holds. Theorem 2 is proved. □. 

Remark 1. It is not possible to estimate the rate of convergence in (8) without 
making additional assumptions on y or on /. In ^ one can find examples illustrating 
similar statements concerning various methods for solving ill-posed problems. 
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